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Abstract 

The aim of this paper is to show some examples of matrix- valued orthogonal functions 
on the real line which are simultaneously eigenfunctions of a second-order differential 
operator of Schrodinger type and an integral operator of Fourier type. As a consequence 
we derive integral equations of these functions as well as other useful structural formulas. 
Some of these functions are plotted to show the relationship with the Hermite or wave 
functions. 



1 Introduction 

In this paper we will show examples ol NxN matrix- valued orthogonal functions ($n)n which 
are simultaneously eigenfunctions of a second-order differential operator of Schrodinger type, 
i.e. 

(1.1) {^nV){x) = ^':,{x)-<S>n{x)V{x)=Tn<S>nix), X £ R, 

where the matrix-valued potential V{x) (independent of n) is a diagonal quadratic matrix 
polynomial (but not a scalar multiple of the identity), and an integral operator of Fourier 
type, i.e. 

1 

(1.2) ($„X)(x) = ^/ <^nit)K{x,t)dt = An<^nix), X G M, 

where the matrix-valued kernel K{x, t) is also diagonal (but not a scalar multiple of the 
identity) of the form K{x,t) = e^^^K for some diagonal matrix K. We will also show that 
it is possible to construct suitable families ($n)n such that the corresponding eigenvalues r„ 
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and A„ are diagonal. Hence both differential and integral operators commute in the space of 
all matrix-valued functions spanned by (<!>„)„. 

Observe that the operators D and I appear on the right. This denotes that the potential 
coefficient and the kernel are multiplied on the right, respectively, while the eigenvalues appear 
on the left. With this configuration it is straightforward to derive that both differential and 
integral operators commute in the space of all matrix-valued functions spanned by ($n)n, i-e. 
for any F = ZZo ^n^n, m > 0, C„ G C^x^, 

FIV = AnFV = KTnF = TnKF = TnFZ = FVI, 

since the eigenvalues are diagonal matrices. In particular we will see in Section [6] that this 
commutativity property holds in C^^^) (see Section[2]for definitions) for the examples 

we consider in this paper. 

In the scalar situation the examples of simultaneously eigenfunctions of an integral and 
a differential operator are very special and usually are related to many areas of engineering, 
physics, chemistry and mathematics. Prominent examples are the Hermite or wave functions 
which satisfy the Schrodinger equation (11.11) with V{x) = and also are eigenfunctions of 
the Fourier transform (see \i?>\ l28l I30j). The Hermite functions play an important role in 
quantum mechanics to model the one dimensional quantum harmonic oscillator (see [23]), or 
in chemistry in vibration of molecules (see [24^ I31j). among other applications. 

Also the prolate spheroidal wave functions in signal processing satisfy certain second-order 
differential equation and are eigenfunctions of the so-called sine kernel. This remarkable 
fact was discovered in a series of papers by the Bell Labs group in the early 1960s, see 
[25l[l9l[20ll26l[27|. See also ^Tj. 

The commuting property and other similar properties are very closely related to the 
classical orthogonal polynomials of Hermite, Laguerre, Jacobi and Bessel. For instance, in 
[15j . the author produces naturally appearing global operators that happen to commute with 
properly chosen local operators. 

In the last few years many examples of matrix-valued orthogonal polynomials {Pn)n have 
appeared satisfying second-order differential equations of Sturm-Liouville type, i.e. 

P'^{X)F2{X) + P'n{x)Fi{x) + Pn{x)Fo{x) = r„P„(x), 

where F2, Fi and Fq are matrix polynomials (which do not depend on n) of degrees less than 
or equal to 2, 1 and 0, respectively, and r„ are Hermitian matrices (if the family {Pn)n is 
orthonormal) . These examples are the matrix analogue of the classical families of Hermite, 
Laguerre and Jacobi polynomials. Since the classical families are intimately related with 
integral equations, it is natural to expect that these new families of matrix-valued orthogonal 
polynomials are related to matrix-valued integral equations. The main goal of this paper 
is to show some examples of simultaneously eigenfunctions of an integral and a differential 
operator in the matrix case. 

For that purpose we will consider families of matrix-valued orthogonal functions (con- 
structed from the matrix-valued orthogonal polynomials and a weight matrix) such that 
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they are orthogonal in the corresponding function spaces (more on this in Section [2]) . This 
treatment will allow us to obtain, for the examples we study in this paper, matrix-valued 
Schrodinger operators of the form (jl.ip and integral equations of the form (jl.2p . 

Schrodinger operators with Hermitian matrix-valued potentials (usually with some re- 
strictions on the potential) are not new in the literature. They have been considered in 
[31 HI nil [22] to study many analytical properties of these operators. Also matrix-valued 
functions are considered for some problems in wavelet theory to model matrix-valued signals, 
e.g. video images (see [291 [32]). 

To give an idea of the results contained in this paper let us display here one of the examples 
we study in Section [H Consider a family of matrix- valued functions {^n)n defined in ()4.3p . 
Then they satisfy the following matrix-valued Schrodinger equation 

<(x) - ^n{x){x^I + 2J) + ((2n + 1)/ + 2J)^n{x) =0, X G M, 

where J is the diagonal matrix defined in ()2.4p . This differential equation (for the monic 
family) was already derived in Section 6.2 of [11] (with non diagonal eigenvalue). The family 
we consider in (j4.3p will transform the eigenvalue into another convenient diagonal one. 
Then the family of matrix-valued orthogonal functions (^n)n satisfies the following integral 
equation 




$„(t)e^^*e*5'^dt = (i)"e^5-^$„(x), x G M, 



where e^2-^ is the diagonal matrix (j2.8p . 

Therefore the kernel and the eigenvalue in (|1.2p are K{x, t) = e*^'*e*2 and A„ = (i)"'e' 2 
respectively. Observe that they are diagonal, but not necessarily scalars multiple of the 
identity. This result will have important consequences. In particular we will derive some 
structural formulas and integral equations of the corresponding matrix-valued orthogonal 
polynomials. 

The paper is organized as follows: in Section [2] we give some preliminaries. In Section 
[3] we derive the matrix- valued differential equation satisfied by the matrix- valued functions 
constructed from the matrix-valued orthogonal polynomials and a weight matrix. We will 
focus on the case when the leading coefficient F2 is of the form F2{x) = f2{x)I, where /2(x) 
is a scalar polynomial of degree less than or equal to 2 (the case considered in [ID]). This 
result holds for any family supported in any interval [a, 6], —00 < a < b < 00. 

In Sections HI and [5] we study in full detail two examples supported in M, which previously 
appeared in [lOj and derive differential and integral equations for the corresponding families 
of matrix-valued orthogonal functions. As a consequence we will derive real (and complex) 
integral equations for the families of matrix-valued orthogonal polynomials, among other 
structural formulas of general size N x N. We also study, for both examples, the special case 
oi N = 2. In this case we can write our (normalized) matrix- valued orthogonal functions 
in terms of classical Hermite or wave functions. We will plot some of the entries of ^n^n 
for the first values of n, as well as derive some integrals involving these functions. Finally, 
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Section [6] addresses the issue of the completeness of these functions, gives a summary of the 
results of the paper and the challenges that lie ahead. 



2 Preliminaries 



An N X N matrix-valued function F is a matrix with real- valued functions as entries of the 
form 



F{x) 



(hi{x) 

/2l(x) 



/12(X) 
/22(X) 



flN{x)\ 
f2N{x] 



X G [a, b]. 



\fNl{x) fN2{x) ••• fNN{x)j 

We will say that F G L'^{[a,hlC^''^) if 

/ F{x)F* {x)dx <oo, 

J a 

where F* denotes the Hermitian conjugate of the matrix F and — oo<a,6<oo. In the 
above definition we mean that the integral is finite entry by entry. It is easy to see that 
F G L^([a, 6], C^^^) if and only if every entry Fij G LF'{a,h). This induces a matrix-valued 
inner product for any two matrix- valued functions F,G^ L^([a, 6], C-^^-^), denoted by 



(2.1) 



{F,G) = I F{x)G*{x)dx. 

J a 



This is not an inner product in the common sense, but it has properties similar to the usual 
scalar inner products. It is also possible to define a scalar product between two matrix-valued 
functions (see |6|), given by 

{F,G) = Tr{{F,G)). 

Therefore, L'^{[a,b],C^''^) with the norm ||F|| = Tr((F,F))^/^ is a Hilbert space and (IZTI) 
is the inner product (in fact, the set of equivalence classes F ~ G if ||F — G|| = 0). Then we 
can module Fourier expansions of orthonormal systems in L^([a, 6], C^^^) (see pp. 7-8 in 

In a similar way, we can define the weighted spaces L^([a, 6], C^^^) of ah matrix-valued 
functions in one variable with the inner product 



(2.2) 



{F,G)w 



F{x)W{x)G*{x)dx, 



where dW is a weight matrix with a smooth density W with respect to the Lebesgue measure, 
satisfying (1) W{B) is positive semidefinite for any Borel set G M, (2) W has finite moments 
of every order and (3) {P, P)w is nonsingular if the leading coefficient of a matrix polynomial 
P is nonsingular. Condition (3) above is necessary and sufficient to guarantee the existence of 
a sequence of matrix polynomials orthogonal with respect to (j2.2p of degree n with nonsingular 
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leading coefficient, and it is fulfilled, in particular, when W is positive definite at an interval 
of the real line. We will assume, for simplicity, that there exist a real number c G (a, b) such 
that W{c) = I. 

Now we introduce two matrices that are going to play a very important role in the rest 
of the paper. The first one is the N x N nilpotent matrix (i.e. A^=0) 

N 

(2.3) A = Y,^jEj,j+i, i^jeW, 

i=i 

where Ej /^ is a matrix with 1 at entry {j, k) and elsewhere, while the second is the diagonal 
matrix 

N 

(2.4) J = ^{N-j)E,j. 

A and J satisfy the algebraic relation 

adA J = [A, J] = -A, 
where ad'^. A; > 0, denote the usual adjoint operators defined by ad^^^ = Y, 

(2.5) adxY = XY-YX and ad^+^Y = adx(ad^y), k > 1. 
As a consequence 

(2.6) ad^fcJ = -A;^^ A; = 1, . . . , iV - 1, 
as it is easy to see by induction. 

We will now introduce real or complex valued functions with matrix argument. The 
Taylor series of an infinitely differentiable / in a neighborhood of x = is 

oo j 

/(x) = 5]/(^)(0)^. 

Whenever we write f{X), for any N x N matrix X, we will mean the following matrix 

fix) = Y,f^'\o)—. 

j=o ■'■ 

It is clear that f{X)g[X) = g(X)f{X) for any two real or complex valued functions f,g. 

In particular we will be interested in the evaluation of f{A) or /(J), where A and J are 
defined in (12. 3p and (12. 4p respectively, for certain real or complex valued functions /. Observe 
first that, since A is nilpotent, then 

fiA) = /(^^(O)-, 
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i.e. f{A) is a finite sum of linear combinations of the powers of A. It is clear that all real or 
complex valued algebraic or differential manipulations for functions f{x) can be applied to 
the matrix f{A). For instance, when we write (/ + A)~^, it will denote the matrix 

N-l 

(i+Ar' = j2i-'^y^'- 

j=0 

On the contrary, the diagonal matrix /(J) will not be a finite sum of linear combinations of 
the powers of J. That is the case of the following diagonal matrices. For f{x) = e^^^^ , k E Z,, 
the matrix /(J) is defined by 

(2.7) e'^''-^ = (i)'^^ 



N 



i=i 



where i is the imaginary unit. These matrices will play a very important role in the integral 
equations that we will study in Sections [4] and [5l Observe that for /c = in (12. 7p we have the 
identity matrix /. For A; = 1 we have the complex diagonal matrix 



N 



(2i 



i- I 



V 



\ 



V 



and for k = 2 we have 



N 



(2.9) 



AttJ 



/(-i)^-i 

V 



V 



which is real and satisfies e™ e™ = I. For A; = 3 we have the inverse of (12. Sh . while for A; > 4 
all matrices will reduce to one of the matrices mentioned above. These diagonal matrices are 
going to play the same role that the imaginary unit i plays in the scalar situation. 

Observe that in (j2.7p we are taking the principal value of log. However it is clear that if 
we take any other value we will have that e*("2+2™'^)'^'^ = e^a'^*^ for any m G Z, so it will not 
give any additional information. 

Similarly we can define the following diagonal (and singular) matrices 



sin(^A;j)=l(e*P^ 
2 / 2i 



e 2 



N 



Y^sm^-k{N-j)]Ej,„ 



/c G Z, 



cos(^^k?j =l(e*f'^^ + e-^f^^) = ^cos(|fc(iV-j))i?,,„ keZ. 
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In the special case of A; = 1 the exphcit expression of sin (f J) and cos (f are the diagonal 
real matrices 



/sin(f(iV_l)) 



(2.10) 



and 



(2.11) 



sin ( —J 
2 



/ 



/cos(|(iV-l)) 



cos ( 



V 



1 / 



This is the only significant case, since for k = 2, sin(7rJ) = and cos(7rJ) = e^'^'^ and for 
k = 3, sin (^J) = - sin (f j) and cos (^J) = cos (f j). 
We remark the following relations 



(2.12) 



e'2-^cos ( — J 



AttJ 



-il + e 



iTTj\ 



e'z'^ sin i — J 



(e^^-^ - /), 



2i 



e'"" cos I "^"^ j ~ (2^'^ 



TT 



e'^i-^ COS i - J 



1 



(/ + 



e""" sin I — J 



e 2 sm —J 



sin ( |j ), 



1 

Yi 



{e'^-^ -I). 



We will use these matrices to derive real integral equations of matrix-valued orthogonal 
polynomials in Section HI 



Now we will prove the following lemma, which gives commutativity relations between the 
powers of A and e^^^^ : 

Lemma 2.1. The matrices A and e^^^'^ , defined in (12. 3p and (12. 7|] respectively, satisfy the 
following algebraic relation 



(2.13) 



giffcj^m ^ (•)ftm^mg*-/cj^ m = 1, 2, . . . , - 1, A; G Z. 



\km Am A^kJ 



Proof: It is enough to prove this for m = 1. For the rest of values of m we apply 
recursively the formula for m = 1. As a consequence of the definitions of A and e^^^'' in ()2.3p 
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and (|2.7p respectively, we have 

N-1 X / AT N N-1 



i=i ^ ^ «=i ^ j=i 

N-1 ^ N X ^TV-l ^ 

i=i ^i=i ^ ^ «=i ^ 

using that Ej^^Eh^m = -^i,m for /c = /i or otherwise. □ 

3 Matrix-valued orthogonal functions 

As we mentioned in the Introduction, in the last few years many families of N x N matrix- 
valued orthogonal polynomials {Pn)n have been found along with their orthogonality measure 
W satisfying second-order differential equations of the form 

(3.1) Pn{x)F2ix) + P;(x)Fi(x) + Pnix)Foix) = r„P„(x), 

where F2,Fi and Fq are matrix polynomials (which do not depend on n) of degrees less 
than or equal to 2, 1 and 0, respectively, and r„ are Hermitian matrices (if the family 
{Pn)n is orthonormal) . These families are natural orthogonal systems in the weighted spaces 
L2^([a,6],C^x^). 

Typically the weight matrices with this property can be factorized in the form 

W{x) = p{x)T{x)T*{x), 

where p is a scalar function (Hermite, Laguerre or Jacobi weight) and T is a matrix-valued 
function which satisfies a first order differential equation with initial conditions of the form 

(3.2) T'{x) = G{x)T{x), T{c) = I, 

for some c G {a,b). The differential coefficient G in (j3.2p is connected with the coefficients 
F2 , Fi of the differential operator (|3.ip and p (see the proof below of Theorem 13.11 for more 
details or |10j). Many examples have been found in the last years by solving a set of three 
symmetry differential equations that are equivalent to the second-order differential operator 
(|3.1j) being self-adjoint with respect to the inner product (j2.2p (see [101 ESI [HI El E]). For a 
different approach, using matrix- valued spherical functions, see |17t I23j. 

One possibility for finding orthogonal systems in the space L^([a, 6], C^^^) is considering, 
for each family of matrix- valued orthogonal polynomials {Pn)n the following family 

(3.3) $„(x) = p^/\x)Pn{x)T{x), n>0. 



Then {^n)n will be orthogonal with respect to the identity matrix /. Considering (13. ip 
and (13. 2p one verifies that the family ($n)ra also satisfies very special differential equations. 
We will restrict ourselves to the special case where the differential coefficient F2{x) in ()3.ip 
is scalar, i.e. of the form F2{x) = f2{x)I, for some real polynomial /2(x) of degree less than 
or equal to 2. This is the case considered in |10] . Therefore we have the following 
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Theorem 3.1. The family {^n)n defined by (j3.3p . satisfies the following second-order differ- 
ential equation 
(34) 

/.(x)<(x) + - (i ( (AW^ + [f^!^) ') / + x(x)) = r„a.„(x), 

where xi^) ^he Hermitian matrix-valued function 

X{x) = T-\x) {j2{x)G\x) + f2{x)G\x) + ilM^KG{x) - Fo)t(x), 



G is the coefficient in (|3.2p and Fq and r„ are the independent coefficient and the eigenvalue, 
respectively, of the differential equation ()3.ip . 



Proof: Differentiating twice (13. 3p one gets expressions of p^/'^P^T and p^^'^P^T in terms 
of the derivatives of ^n- Multiplying the differential equation (j3.ip on the right by p^l'^T and 
substituting the formulas mentioned above we get (13. 4p after using Fi = 2/2G + ^•^^^ (see 
formula (4.5) in [lOj). The matrix- valued function x is Hermitian as a consequence of the 
third symmetry equation (see formula (4.12) in [10] for details). □ 

Observe that the coefficients of the differential operator ()3.4p are extremely simplified and 
are all scalar functions except for the independent coefficient or potential, which is Hermitian. 
In fact, in all the examples in the literature until now, the matrix x is a diagonal matrix (see 
comment in pp. 93 of [9j). 

The expression of the differential equation (j3.4p for each one the classical weights of 
Hermite, Laguerre and Jacobi is 

1. For p{x) = e~^^ and f2{x) = 1 

(3.5) <(x) - $„(:e)((x2 - 1)/ + x{x)) = TAix). 

2. For p{x) = x°'e~^ , a > —1, and f2{x) = x 

(3.6) X<(X) + Kix) - <fn(x)((j - ^ + + = ^nM^)- 

3. For p{x) = (1 - x)"{l + xf,a, (3 > -1, and /2(x) = 1 - 
(3.7) 



l-x2X(x)-2xcI>;(x)-^.„(x)((2(fe) + 2(fe) - ("+^)(^+^+^) ^/ + x(x)) = TA{a 



The explicit expression of the Hermitian matrix- valued function x in each case can be found 
in Lemma 2.2 of p] for ([33]) and 1^3^ and in formula (2.5) of for (fHTTp . 

Observe that the equation (13. 5p may be regarded of as a matrix-valued version of the 
one-dimensional Schrddinger equation. The equation (|3.6p is related with a matrix-valued 
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version of the confluent hypergeometric equation, while the equation (j3.7p is related with a 
matrix-valued version of the differential equation for spheroidal wave functions (for special 
values of the parameters a and /3). These differential equations in the scalar case are very 
important with diverse applications to physics, engineering and mathematical analysis itself 
(see [a 1211 [2]). 

In this paper we will focus only on families satisfying matrix-valued Schrodinger differ- 
ential equations as in (j3.5p . i.e. supported in the real line. For each of the families that we 
consider, we will prove that there exists an integral operator of the form (jl.2p with kernel 
K{x, t) having the matrix-valued functions {^n)n as eigenfunctions. We will concentrate on 
Examples 5.1 and 5.2 of |10]. 

We will always denote by ($n)n the matrix-valued functions for all the examples, while 
there is no confusion about which family we are using from the context. In the case we are 
referring to different examples at the same time we will denote them by {^n.k)ni k = 1,2, 
where {^n,i)n is the example in Section [J] and ($n,2)n is the example in Section [5l 

4 The first example 

Let W be the following weight matrix 

(4.1) Wix) = e-^'e^^e^*^', xGR, 

where A is the N x N nilpotent matrix (j2.3p . 

This example was considered for the first time in [lOj. The family of monic orthogonal 
polynomials {Pn)n satisfies a second-order differential equation as in (|3.ip with 

F2{x) = I, Fi{x) = -2x1 + 2A, Fo{x)=A^-2J, = -2nl + A^ - 2J, 

where J is the N x N diagonal matrix defined in (j2.4p . 

Now we construct a family of polynomials of the form Pn{x) = LnPn{x) where the leading 
coefficient L„ is chosen such that the eigenvalue r„ transforms into one diagonal. A possible 
choice for this example is 

K = e-^''\ 

It is straightforward to see, using the formula e^He~^ = X]j>o ^dx(^)/.?' (^^^ (|2.5p for 
definitions) and ([Sj]) for k = 2, that L^T^L"^ = -2nl - 2 J. 
Hence, denoting 

(4.2) P4x) = e-^'/^P„(x), 
the family 

(4.3) $„(x) = e-^"'^Pn{x)e''^ 
satisfies a more convenient differential equation 

(4.4) <(x) - ^n{x){x'^I + 2J) + ((2n + 1)/ + 2J)$„(x) = 0, 
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as a consequence of Theorem 13.11 (see (|3.5p ). Observe that in this case we have xi^) = 2 J 
(see Theorem 5.1 of [lO]). We remark that this differential equation is independent of the 
matrix A. Note that a similar differential equation (using the monic family) was derived in 
Section 6.2 of [H]. 

In order to prove the main result in this section (Theorem 14. 1 1 below) we need the following 

Lemma 4.1. The following formula holds 
1 

(4.5) / e-*'/2e-AV4eAig«igif j^^ ^ jg-.V2e-AV4gAx^ 

V27r J-oa 

where e*2"^ is the diagonal matrix (12. Sp . 



Proof: Expanding e^* on the left hand side of the formula (|4.5|) . using formula (|2.13|) 
for k = 1 and e~'^^/^e*2'' = e*2"^e'^^/'^ (as a consequence of (|2.13p ). and denoting Hn{x) = 
(— l)"e^ /^(e~^ f^iiQ monic Hermite polynomials, we obtain 



1= re-*'/2e-^'/V*e-*e*i^di=^y e-^'/^^^e^f-^ H e-'^'^^'^'dt 

\/2^t:;^ J-oo j! 



= ^ e^l-'e^'/4^(-i)^e-^'/2(i)^-^,(x) 
i=o -^^ 



since the monic Hermite functions Hn{x)e ^^/^ are eigenfunctions of the Fourier transform 
with eigenvalue (i)" and the generating function for the monic HermitG polynomicils (^Hfi^n 
is given by Hjix^j = e^*-*V2. □ 

Theorem 4.1. The family of matrix-valued orthogonal functions {^n)n defined in ()4.3p sat- 
isfies the following integral equation 

(4.6) / $„(t)e*^*e^i-^dt = (i)"e^f ^$„(x), x e M, 

where e*2-^ is ^/le diagonal matrix (|2.8p . 
Proof: Denote 

(4.7) *„(x) = ^ / ^nit)e'^'e'^-^dt. 
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By integration by parts using ^e*^* = zte*^* and ^e*^* = ixe"* we get that '^nix) satisfies 
the same second-order differential equation as ^'^(x), i.e. ()4.4p . Therefore ^nix) can be 
written as "^n{x) = C„<l>„(x), n > 0, for some sequence of nonsingular diagonal matrices Cn- 
This follows if we expand "^nix) = YlT=o^'n',k^k{x)- Since ^nix) satisfies (|4.4p and ($n)n 
is a system of linearly independent matrix-valued functions this is equivalent to say that 
(n — k)Cn,k + JCn,k — Cn,kJ = for all n,k > 0. But now it is straightforward to conclude 
that if n ^ k then C„^fc = and if n = A; then Cn,n = Cn must be a diagonal matrix, since J 
has simple spectrum. 

Now we will use Lemma [4. 11 which is exactly the case n = in (j4.7p . After differentiating 
the right hand side of the formula ()4.5p n times with respect to x one obtains 

(4.8) e^f ^e-^ /^-^(e-^ /^e^^) = ^ ^ % / e'' /H^e^'e'^'e'^'dt. 



Using the Leibniz's formula, the left hand side of the formula (14. Sp can be also written as 
(4.9) e*f ■^e-^'/^£^(e-^'/2g^x) ^ e*f ^e-^'/^((-l)"x"/ + • • • )e.--^l''e^\ 

Writing <l*n(x) = e~^^/^P„(x)e'^^ = e~^^/^e~'^^/^(x"/ )e"^^, by linearity, and using the 
right hand side of (14. 8|) and (14. 9|) one obtains 

^'„(x) = ^ / «>„(t)e^^*e^5-^di = / {fl + . . . )e^*e^^*e^i-^di 

v27r J-oo v2vr y_oo 

= (i)"e^f ■^e-^'/^(x"/ + • • • )e-^'/2gAx_ 

Since we already know that ^n{x) = Cn^n{x)., equating the leading degree of x in the formula 
above we obtain C„ = {i)^e^2-^ . □ 



Observe that the integral equation (14. 6p is also independent of the matrix A. 
The integral equation will have important consequences: 

Corollary 4.1. The family of matrix-valued orthogonal functions ($n)ra satisfies the following 
symmetry condition 

(4.10) $„(x) = {-ire'^^<^>ni-x)e'''\ 

where e™'^ is the real diagonal matrix (12. 9p . Consequently, the family of matrix-valued or- 
thogonal polynomials Pn{x) = e~^ I Pn{x\n > 0, satisfies the same symmetry condition, 
i.e. 

(4.11) Pn{x) = {-ire'^'Pn{-x)e'^'. 
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Proof: From (j4.6p multiplying on the left by the eigenvalue (i)"e*2-^ and substituting 
again by the same formula we get 

-| poo 1 poo poo 

V27r y_oo ^vr y„oQ 

/oo 
-oo 

using standard Fourier analysis, where S is the standard Dirac delta function. Therefore 
(j4.10p holds. The formula (j4.1ip holds from the observation that e~^^e*'^'^e~'^^ = e^'^'^ , as a 
consequence of (j2.13p . □ 

Corollary 4.2. The family of matrix-valued orthogonal polynomials Pn{x) = e~^^^^Pn{x), 
n > 0, with respect to the weight matrix (14. ip satisfies the following integral equation 



(4.12) e-^'/^Pn{x)e^^ = ^-^e-'^' H Pnit)e~'' /^e'^'e^'e'^' dt, 

V27r J~oo 

where e*2'^ is the diagonal matrix ()2.8p . 

Proof: Immediate from (|4.6p using (|4.3p . □ 



We remark here that the family {Pn)n is a solution of the integral equation ()1.2p with 
kernel = e(^+^*)'e^*e*f ^e"^^. 

Observe that (j4.12p is an integral equation which depends on complex values. In the scalar 
case it is well known that the Hermite polynomials {Hn)n satisfy real integral equations in 
terms of the kernels cos(xi) and sin(xi) (see for instance [21j). This is possible since the 
Hermite polynomials satisfy the symmetry condition Hn{x) = {—l)'"'Hn{—x). In our case, 
we have a different symmetry condition (j4.1ip . so it will not follow the same lines as in the 
scalar case. However, it is possible to derive real integral equations for {Pn)n- 

Corollary 4.3. The family of matrix-valued orthogonal polynomials Pn{x) = e~^^^'^Pn{x), 
n > 0, with respect to the weight matrix (|4.ip satisfies the following real integral equations 

(4.13) e-^'/2(e*"-^ ± /)P„(x)e^^C± = ^ ' C± / e-*'/2fc„(a;, t)P„(t)e^*dt(e*"-^ ± /), 

\/2tT J-oo 



and 
(4.14) 



e-V\e'-'±I)Pn{x)e^''C^ = ^^^t^C± H e-''/^kn+i{x,t)Pr.{t)e^'dt{e'^'' T I), 

v27r J-oo 

where e^'^'^ is defined in (j2.9p . C_ = sin (f J), C+ = cos (f and are both defined in ()2.10p 
and ()2.1ip respectively, [-J denotes the floor function and 



kn{x, t) 



cos(xt), if n is even, 
sm{xt), if n is odd. 
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Proof: From the symmetry condition (j4.1ip and using formula (j2.13p we have that 
Therefore, the evaluation of (j4.12p at x and —x and using the previous formula gives 



Now we multiply on the right and on the left by appropriate matrices such that the elements 
in the big parenthesis of the formula above are equal and the other vanishes for even or odd 
values of n. This will become clear below. These matrices are sin (§ J) and cos (f defined 
in (j2.10p and (j2.1ip respectively. There are four possible combinations of multiplying these 
matrices on the right and on the left. This is why we will get eventually eight formulas, 
considering even and odd values of n. We will show the case when we multiply on the left by 
cos (f ^) and on the right by sin (f J). From the analysis below the rest of formulas can be 
derived in a similar way. 

From the previous formula and using relations (I2.12P one gets 



e -'/2(e--^ + /)P„(x)e^" sin i-j] = cos f- J 



oo 

*V2 ^■r.^(r,.+\((_^\n _ -\\p_{f\^At 



2 J 2iV2^ \2 

/oo 
e-*'/2 sin(xt)((-l)" + l)Pn{t)e'^'dt 
-oo 



e-* /"cos(xt)((-l)" - l)P„(t)e^*dt 



Now it is clear that for even or odd values of n, either the integral with cos(a;t) or sin(xi) 
will vanish and that these integral equations are real. Therefore 

e-"'/2(e^--^+/)P2n+i(a;)e^^sin = cos (^^j) J"" e-'" cos{xt)P2n+i{t)e^'dt{e'^' -I), 

which are formulas (j4.14p for the positive sign and both odd and even values of n. □ 

Observe that (|4.13p and (|4.14p give eight real integral equations of the family {Pn)n 
according to positive or negative sign and odd or even values of n. The reason for eight 
formulas is the structure of the matrices sin (f J) and cos (f J), defined in (|2T0]1 and (fZTT]) 
respectively. These are diagonal singular matrices and whenever we multiply them on the 
right or on the left we only get a description of some of the entries of P„ (which is a full matrix 
in general) . It is exactly this combination of multiplying on the right and on the left by these 
matrices that allows us to have a formula for every entry of P„- In other words, if we multiply 
on the left by cos (f ^) and on the right by both cos (f or sin (f J), we get a description 
of the N, N — 2, N — A, . . . rows of Pn while if we multiply on the left by sin (f J) and on the 
right by both cos (f ^) or sin (f J) , we get a description of the A^ — 1,A^ — 3, iV — 5,... rows 
of Pn, respectively. Therefore all rows of Pn are covered. 
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Remark 4.1. We have found that the family of polynomials {Pn)n defined in ()4.2p (and 
the family we will study in Section defined in (|5.2p ) simplifies considerably many structural 
formulas. For instance, the norms of {Pn)n (and consequently the norms o/ with respect 
to the inner product (j2.ip ) are diagonal. Also the coefficients of the three-term recurrence 
relation are considerably simplified. A detailed study of these and other structural formulas 
for {Pn)n will be discussed in future publications. 



4.1 A detailed study of the case N = 2 

In this section we will study in detail the properly normalized matrix-valued functions ($n)n 
for the special case of = 2. We will show their relationship with the scalar Hermite or 
wave functions, derive some structural formulas and plot graphs of the diagonal entries of 
for some values of n and ui (the only free parameter in the matrix A, see 



Using Theorem 5.1 of JT] one can derive an explicit expression for the family of matrix- 
valued orthogonal polynomials {Pn)n defined in (14. 2p . They can be written as 

p ( ^ _ ]_ f Hn{x) -nuiHn-i{x) \ 

^n[X) - 2„ [^_n^Hn-l{x) l{Hn{x) + nufxHn-l{x)) J ' 

where Hn{x) = {—l)^e^ (e~^ )'"', n > 0, are the classical Hermite polynomials, and (7„) is 
a sequence of real positive numbers defined by 

(4.15) 7n = 1 + 2^?' ^ > 0- 

We can normalize this family since the norms of {Pn)n are diagonal 

IIP ||2 _ fjn+i 

W^nWw- 2" V 1/7. 

Therefore the normalized matrix- valued functions = ||P„||^^<I>„ can be written as 



V'n(x)/^/7^ i^lY^2^V'n+l(x)^ 
-^-^J^^n-lix) ^n{x)l^ 



(4.16) ^nix) 



where i^nix) = , ^ e ^'^^'^Hn{x) are the normalized Hermite or wave functions. Hence 
(4.17) 

^ ^ ^ ^ 7^feV'n(x)V;(x) ^l;l__^{x) + ^{i;l{x)-^l^l_^{x))J ' 

In the expressions above we are using standard properties of Hermite polynomials and wave 
functions like Hn+i{x) = 2xHn{x) — 2nHn-i{x) and i(>'n{x) = \J~^^'4)n+i{x) — ^J^^^l>n-l{x). 
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Observe from the definition of (7n)n in (|4.15p that the diagonal entries of (j4.17p are probabihty 
densities depending on one free parameter vi. 

The effect of the parameter vi can be studied from the exphcit expression of ()4.16p and 
(j4.17p . For instance, as z^i tends to we observe that the diagonal entries of ()4.16p converge 
to two copies of the Hermite functions while the off diagonal entries tend to 0. If ui is large, 
the diagonal entries get small compared with the off diagonal entries (for the first values of 
n). 

In Figure 1 and 2 we have plotted the diagonal entries of $„(a;)<I>* (x), which are prob- 
ability densities in M, for every n = 0, 1, . . . , 5, and for the special choice of ui = 1. The 
numbers on the graphics correspond to the value of n. It is well known that in the scalar 
situation there exist n points (the zeros of the Hermite polynomials) where the probability 
density tpn{x)tp^{x) is exactly 0. In our case our two probability densities never vanish. 




Figure 1: ($„$;)ii, n = 0, . . . , 5, z^i = 1 



Figure 2: {<^nK)22,n = 0, . . . , 5, z^i = 1 



Finally we compute the following formulas 



(4.18) 



x''$„(x)^>;^(x)(ix, n,m>0, k = 1,2. 



In the scalar situation, we have, as a consequence of the three-term recurrence relation among 
Hermite polynomials, that 



(4.19) 



'm,n— 1 



+ 



n + 1 



(4.20) 



{x^)nm = -\/n(n - l)6m,n-2 + {n + l/2)5m,n + -yj {n + l){n + 2)5m,n+2- 
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We can use the formulas in the scalar case to calculate the matrix- valued formulas (j4.18p . 
From ()4.16p we can obtain an explicit expression of <I>n^m- 



^n^Pm I / (n+l)(m+l) , 

V7n + l7m+l 2 '' 



^ U In ym 



m—1 



And using formulas ()4.19p and (|4.20p for the scalar wave functions we obtain 





ixl)r 








and 



1 / n(n- l)7n+i 

2 V 7n-l 



+ 



+ 



771 + 1 

n 




v7n 



7n + l 



27n-l7n+l I A 



n+1 



71+1 

27n7n + 2 



-'m,n+l 



1 / (n+l)(n+2)7n+3 

2 V 7n+l 





1 / (n+l)(n+2)7 

n 

2V TwTi 



+2, 



where (7n)n is the sequence (j4.15p . Observe that for A; = 1 we get an extra term when n = m, 
something that did not happen in the scalar situation. Likewise for k = 2 we get two new 
extra terms. 

In the scalar situation the tridiagonal matrix (x) given in ()4.19p can be viewed as the 
matrix of the homomorphism / i— )• x/ in L^(M) with respect to the basis {ipn)n- In the 
matrix case for this example we have that the block tridiagonal matrix (xl) is the matrix of 
the homomorphism F i— )• xF in L^(M, C^^^) with respect to the basis (^*n)n- In this case, 
(xI) is a 4 diagonal semi-infinite matrix with zeros in the main diagonal of the form 



(x/) 



/O * * 

* 

* 

V 








where a * means a nonzero entry. 
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5 The second example 

Let W be the following weight matrix 

(5.1) = e-"'e^(^+^)"'"%(^+^*)"'^*"', xEM, 

where A is the N x N nilpotent matrix (j2.3p . This weight matrix was considered for the first 
time in Example 5.2 of [10], where the authors use the notation B = A{I + A)~^. Although 
this new notation looks more complicated it will simplify considerably all computations in 
this section. All proofs follow the same lines as in Section HI 

The family of monic orthogonal polynomials {Pn)n satisfies a second-order differential 
equation as in (13. ip with 

F2{x) = I, Fi{x) = 2{2 A{I + A)-^ - I)x, Fq{x) = 2A{I + A^^ - 

Vn = 2{2A{I + Ay^ - I)n + 2A{I + - 4J, 

where J is the N x N diagonal matrix defined in (j2.4p . 

Now we construct a family of polynomials of the form Pn{x) = LnPn{x) where the leading 
coefficient L„ is chosen such that the eigenvalue r„ transforms into one diagonal. Observe 
that the eigenvalue has n dependence off the main diagonal. Therefore the coefficient L„ will 
depend on n. A natural candidate in this case is 

as a consequence the algebraic relation \og{I + A)J — Jlog(/ + ^) = —A{I + A)~^ which can 
be proved using (12. 6p (just expanding in power series). 
Hence, denoting 

(5.2) Pn{x) = [(/ + A)-i/2]2"+ip„(x), 
the family 

(5.3) ^n{x) = e--'/2p„(x)e^(^+^)"'-' 
satisfies a more convenient differential equation 

<(x) - $„(x)(x2/ + 4J) + ((2n + 1)/ + 4J)$„(x) = 0. 

as a consequence of Theorem 13.11 (see (|3.5p ). Observe that in this case we have xix) = 4 J 
(see Theorem 5.2 of [lOj). We remark again that this differential equation is independent of 
the matrix A. Note that a similar differential equation (using the monic family) was derived 
in Section 6.2 of [II]. 

As before, in order to prove the main result in this section, we need the following 
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Lemma 5.1. The following formula holds 

1 poo 

where e^'^'^ is defined in (j2.9p . 



Proof: Although the proof is more elaborated than the one in Lemma [4. 11 it follows the 
same procedure using standard real analysis of power series of functions and the formulas 
gmj^fcj^j + A)-^f = {-!)'' A''[{I + e'^-^ and (/ + A)V2gmJ ^ ^mJ^j _ ^)i/2^ ^-^ich 

hold using Lemma l2.1[ □ 

Therefore we have the following 

Theorem 5.1. The family of matrix-valued orthogonal functions {^n)n defined in (I5.3P sat- 
isfies the following integral equation 



1 f-OO 

V27r J-oo 



where e^'^'^ is defined in ()2.9 



Proof: The only difference with respect to the proof of Theorem 14. II is that the leading 
coefficient of P„ depends on n and that the formula (j4.9p is now 



e 



dx"^ 

= ((-1)"[(/ - A){I + Ay^Y^x""! + ■■■ )e-^V2^Mi+A)-^x^y 

The rest follows the same arguments as in Theorem 14.11 □ 

In a similar way we have the following corollaries, from which we omit the proofs since 
they are exactly the same as in the example in Section [H 

Corollary 5.1. The family of matrix-valued orthogonal functions ($n)n o-'f^d the family of 
matrix-valued orthogonal polynomials {Pn)n satisfy the following symmetric conditions 

(5.4) <D„(X) = {-ir^n{-x), Pn{x) = (-l)"Pn(-x). 

Observe now that the symmetry conditions are exactly the same as the classical Hermite 
polynomials. 

Corollary 5.2. The family of matrix-valued orthogonal polynomials {Pn)n, n > with 
respect to the weight matrix (15. ip . satisfies the following integral equation 

e-^-'/2p„(x)e^(^+^)"'"' = ^Le*'^^ / Pr,{t)e''' /^e'^'e^^'+^^'"'\'^J dt, 

v27r y_oo 
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where e^'^'^ is the diagonal matrix (j2.9p . Moreover, we have the following real integral equa- 
tions 



OO 



g-xV2gi-Jp2n+i (rE)e^(^+^)"'^' = / e-*'/2p2„+i (t)e^(^+^)"*' sin(xt)e^"'^dt. 

V2vr J^oo 

Observe now, because of (|5.4|) . that there are only two real integral equations, unlike the 
example studied in Section [H 

5.1 A detailed study of the case N = 2 



Using Theorem 5.1 of [11] one can derive an explicit expression for the family of matrix- valued 
orthogonal polynomials {Pn)n defined in (|5.2p . They can be written as 

p.^]_( Hn{x) -ui{{n+l/2)Hn{x)+n{n-l)Hn-2{x)) 

where Hn{x) = (— l)"'e^ (e~^ ) ■, n>Q are the classical Hermite polynomials, and (7n)n is 
a sequence of real positive numbers defined by 

2 / ■ 



(5.5) ^" = ^ + f(2j' 

We can normalize this family since the norms of {Pn)n are diagonal 



2 _ ( 7„+2 



1/7. 



Therefore the normalized matrix- valued functions = 'J'n can be written in the 

following form 



(5.6) $„(x) 



MX)/^1^2 ^y^^^|±^^„+2(x) 



where ijjnix) = e ^^/^ff„(x) are the normalized Hermite functions. Therefore 

(5.7) 

ln{x)K{^) - 



^1+2{X) + 7^,{i^l{x) - i^l+2{x)) - 2^"X+, MX){MX) + 'ixil,'^{x)f 
- VWlV'n(x)(V'n(x) + 2X<(X)) ^'-2(2^) + ^{^l{x) " i^l-2{x)) ^ 



Observe from the definition of (7n)n in (|5.5p that the diagonal entries of ()5.7p are probability 
densities depending on one free parameter vi. The effect of the parameter vi is similar to 
the example studied in Section HI 
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Figure 3: ($„$;)ii, n = 0, 5, z^i = 1/2 



Figure 4: ($„$;)22, n = 0, 5, vi = 1/2 



In Figure 3 and 4 we have plotted the diagonal entries of (x) for n = 0, 1, . . . , 5, 

and for the special choice of vi = 1/2. The numbers on the graphics correspond to the value 
of n. 

Finally we compute the following formulas 

/oo ^ ^ 

x^^n{xWm{x)dx, n, m > 0, = 1, 2, 
-oo 

in a similar way as in the previous example using formulas (I4.19P and (I4.20p . 
From (j5.6p we can obtain an explicit expression of $„<5^: 



V'n(3^)V'm(x) / (n+l)(n+2)(m+l)(m+2) . x , 

V 7n + 27 m + 2 ^ ~V 7n + 27m + 2 " + 2 (.X J ^^+2 (.^J 

'^^^^±^^V'™+2(X)V„(X) - ^^^i^m{x)i^n-2{x) 



y v'7n7m + 2 \ 2 



Therefore we obtain 



"Tri+2 \ / / ("+l)7n + 3 / n + 1 

f^n - I V ^^"+1 , W 1 + V 27„+2 ly 27„+i7„+ 2 1 e 



27n7n + l y "^-yn / \ V 27 



Tn + 1 
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and 



1 / n(n-l)7„+2 
2V 7n 
i^i -^n(n-l) 

7n 







6„ 



+ 



1 / (n+l)(n+2)7n+4 

2 V 7n+2 



n(n-l)7„_2 

7n 

7n + 2 



-2 + 



n + 



2 7n + 2 

i^i(2n+l) 

2V7n7n+2 



vi(2n+l) N 

2V7n7n+2 

3 I J. 

2 7n/ 



n 







(n+l)(re+2)7 

7n + 2 



5m,n+2 ) 



where (7n)n is the sequence (jS.Sp . 

Now the block tridiagonal matrix (x/) is the matrix of the homomorphism F 1— )■ xF in 
L^(R, C^^^) with respect to the basis {^n)n- In this case, {xl) is a 7 diagonal semi-infinite 
matrix with zeros in the first three main diagonals of the form 

\ 



{xl) 








★ 









































* 


































\ 



where a :*r means a nonzero entry. 



6 Concluding remarks 

In this paper we have shown the first examples of integral operators having families of matrix- 
valued orthogonal functions as eigenfunctions, which at the same time are eigenfunctions of 
a second-order differential operator of Schrodinger type. We have focused on two examples 
supported in the real line appeared in [1^. Certainly other examples supported in the real 
line have been found as well. For instance, the example in Theorem 1.1 in [9j and the 
example included in the Appendix A.l in [9]. These examples are generalizations of the 
example studied in Section HI We have found that these examples also satisfy similar integral 
equations of the form (11. 2p , along with similar structural formulas and integral equations for 
the matrix-valued orthogonal polynomials. 

These integral operators are slight modifications of the usual Fourier transform, but 
not necessarily the scalar Fourier transform. Instead, we can perform Fourier analysis for 
matrix- valued functions using the families of matrix-valued orthogonal functions introduced 
in Sections [4] and [5l This is different from the usual Fourier analysis where the kernel is 
K{x^ t) = e*^*I, and the transform is applied entry by entry. For instance, for the examples 
^n,k, k = 1,2, introduced in Sections H] and [5] respectively we can recover any matrix-valued 
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function F G L^(M, C^^^) that can be written in the form 

oo 

= ^Cn,k^n,k{x), Cn,k = {F,<^n,k), k = 1,2, 
n=0 

with the inner product (•, •) defined in ()2.ip . This representation is in fact always possible for 
any F G L2(M,C^x^), since our families of eigenfunctions ($n,fc)n; k = 1,2, are complete in 
L2(M,C^'<^). We omit the details, but this is a consequence of having our weight matrices 
(14. ip and (15. ip as nonsingular matrices that die out exponentially at infinity to the zero 
matrix, that is 

Ty(x) = 0(e~"'^') for some a>0 as \x\ ^ oo. 

Hence it is possible to extend to the matrix case all the results necessary to proof the 
completeness of Hermite (and Laguerre) polynomials (see, for instance, the Theorem appeared 
in [18] or Section 6.5 of [Ij) and use Fourier expansions (see pp. 7-8 in [6j) to conclude that 
our families {^n,k)n, k = 1,2, are complete in L^(M, C^^^). 

Therefore, define the following two integral transforms of F 

{FFk){x) = -= / F{t)e'-'e'^^^'dt, k = l,2, 

where e^^^'^ is defined in ()2.7p and Fk acting on the right means that the kernel is multiplied 
on the right. The case /c = is the usual Fourier transform entry by entry but if we multiply 
on the right by e*^^'' we can recover F in a nicer way using the eigenfunctions ($n,fc)n) 
k = 1,2. If we define the inverse of the previous two integral transforms as 

1 f °° 

(FJ--i)(x) = ^ / F{t)e-'^'e-'^''^dt, k = l,2, 

it is easy to see that {<l>n,kJ='k){x) = (i)"e*i''-^^>„,fc(x) and ($„,fc-^^"^)(x) = ( -i)"e ^^^-^^„,fc(3;), 
using the symmetry condition ^n,k{^) = {-l)'^e^'^'^-^ ^n,k{-x)e^'^'''^ (see (!4.10p and (|5.4p ). 
Therefore 

F(x) = {{FFk)F^'){x). 

That means that Fourier analysis of matrix- valued functions in L^(M, C^^^) can be 
studied, at least, in two different ways according to the Fourier type transform that we 
consider and its corresponding eigenfunctions, not only applying Fourier transform entry by 
entry. 
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